Energy-resolved inelastic electron scattering off a magnetic impurity 
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We study inelastic scattering of energetic electrons off a Kondo impurity. If the energy E of the 
incoming electron (measured from the Fermi level) exceeds significantly the Kondo temperature Tk, 
then the differential inelastic cross-section o~(E,uS), i.e., the cross-section characterizing scattering 
of an electron with a given energy transfer u, is well-defined. We show that a{E,ui) factorizes 
into two parts. The independence of a(E,u) is logarithmically weak and is due to the Kondo 
renormalization of the effective coupling. We are able to relate the cj-dependence to the spin-spin 
correlation function of the magnetic impurity. Using this relation, we demonstrate that in the 
absence of magnetic field the dynamics of the impurity spin causes the electron scattering to be 
inelastic at any temperature. At temperatures T low compared to the Kondo temperature Tk, the 
cross-section is strongly asymmetric in u and has a well-pronounced maximum at hu> ~ Tk- At 
T S> Tk, the dependence a vs. oj has a maximum at u — 0; the width of the maximum exceeds Tk /h 
and is determined by the Korringa relaxation time of the magnetic impurity. Quenching of the spin 
dynamics by an applied magnetic field results in a finite elastic component of the electron scattering 
cross-section. The differential scattering cross-section may be extracted from the measurements of 
relaxation of hot electrons injected in conductors containing localized spins. 

PACS numbers: 72.15.Qm, 75.20.Hr, 73.23.-b 
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I. INTRODUCTION 

Scattering of an electron off a magnetic impurity em- 
bedded in a conductor is known to be anomalously 
strongi. The origin of the anomaly is rooted in the de- 
generacy of the localized spin states. This degeneracy, 
being removed by a weak exchange interaction with the 
itinerant electrons in a metal, gives rise to the strong 
scattering of electrons with low energy - the Kondo ef- 
fect. Perturbation theory in the exchange interaction 
constant J is singular. The second-order contribution in 
J to the scattering amplitude diverges logarithmically 2 
if the electron energy E (measured from the Fermi level) 
and temperature T are approaching zero. It is important 
to notice that the logarithmically divergent contribution 
to the amplitude corresponds to an elastic process. In- 
deed, this contribution comes from the change of state of 
one electron; states of all other itinerant electrons are the 
same in the beginning and end of the scattering process. 
Therefore, the energies of the electron before and after 
the scattering is unchanged. 

The divergence noticed by Kondo is not unique to the 
second order of the perturbation theory. Its higher or- 
ders (n > 2) also contain divergent terms of the type 
jn l n rl_1 (£)/ £ ) ] where e — max(E, T), and D is some ul- 
traviolet energy cut-off, whose value depends on the spe- 
cific model; D 3> s. These leading logarithmic terms may 
be summed up by diagrammatic method^ or by means of 
the "Poor man's scaling" 4 renormalization group (RG), 



yielding for the scattering amplitude 



S->k',a'S' 



Hs/Tk) 
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where s CT .' and S S)S ' are the spin operators of the con- 
duction electrons and the impurity, respectively. The so- 
called Kondo temperature is given in terms of the cut-off, 
D, and the exchange interaction, J, as Tk = De~ 1 '^ Ju ', 
where v is the density of states. Like the lowest-order 
perturbation theory result, the leading-logarithmic ap- 
proximation Eq. JQ| corresponds to purely elastic elec- 
tron scattering. 

The leading-logarithmic approximation is adequate at 
e ^> Tk, but it fails at low temperatures. A convenient 
phenomenological description of the low-energy behavior 
of a single-channel Kondo model is given by Nozieres' ef- 
fective Fermi liquid theory. In this theory, a scattering 
problem can be formulated, too. It is clearly seen&, how- 
ever, that the scattering is not purely elastic at e -C Tk- 
At T — 0, for example, the inelastic contribution to the 
electron scattering cross-section scales as (E/Tk) 2 , and 
becomes comparable to the elastic part at E ~ Tk- 

The Kondo effect is a crossover phenomenon, rather 
than a phase transition. The measurable characteristics, 
such as the contribution to the susceptibility or resistivity 
due to magnetic impurities depend smoothly on temper- 
ature. Similarly, the electron scattering off a magnetic 
impurity, which is deeply inelastic at e ~ Tk, must have 
some inelastic component at any energy E. In this pa- 
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per we investigate in detail the inelastic scattering of a 
high-energy electron off a magnetic impurity. 

A study of the energy-resolved, differential cross- 
section, <j(E,uj), is interesting in its own right, but it 
can, in principle, also be measured, e.g., in a modification 
of the experiments of Pothier et al&. Further motivation 
to study o~(E, uj) beyond perturbation theory comes from 
the recent theoretical work of Zarand et alrL. In Ref. the 
energy dependence of the total scattering cross-section, 
o~tot{E) = J duj <j(E,uj), was addressed. With the help 
of the optical theorem, the total cross-section a to t(E) 
was compared with the elastic part of it. The conclusion 
reached in Ref . [7] regarding the energy domain E Tk 
is striking: at T = the scattering is deeply inelastic; the 
elastic part turns out to be negligibly small. This seem- 
ingly contradicts the leading-logarithmic result for the 
scattering amplitude Eq. |T]). The physical explanation 
of this phenomenon however remained unclear in Ref. 7 
and motivates us to revisit the problem of inelastic scat- 
tering. The dependence of the differential cross-section 
o~(E,u>) on uj, which we consider in this paper, clarifies 
the issue, as we are able to determine the distribution 
of energy losses in the inelastic electron scattering off a 
magnetic impurity. 

The separation of the electron scattering cross-section 
in the Kondo effect into elastic and inelastic parts at 
E » Tjc was not addressed for decades, as it does not 
affect the routinely measured quantity, the resistivity. 
The anomalously fast electron energy relaxation in some 
mesoscopic metallic wires&, which was discovered in the 
last decade, prompted a search for relaxation mechanisms 
driven by impurities with internal degrees of freedom. A 
viable mechanism of energy relaxation was suggested first 
in Ref. |8|, and was associated with the electron-electron 
scattering mediated by exchange interaction of electrons 
with magnetic impurities. The removal of degeneracy 
of the localized spin states by the exchange interaction 
results in an anomaly of the electron-electron scatter- 
ing cross-section at small energy transfers^; the colli- 
sion of two electrons with energies E,E' ^S> Tk leads 
to a re-distribution of the energies between the two par- 
ticles, E, E' — > E — ui, E' + u), and has cross-section 
K(u>, E, E') oc J a /uj 2 in the lowest-order perturbation 
theory-^. The l/u> 2 dependence of K allowed the ex- 
perimental observations^ to be explained qualitatively. 
Later experiments^ performed in a magnetic field suffi- 
cient for the Zeeman splitting of impurity energy levels 
did confirm the origin^ of the inelastic electron-electron 
scattering, and indicated the irrelevance of more exotic 
mechanisms, which assumed a generic non- Fermi liquid 
behavior introduced by impurities 1 ^. 

The existence of energy exchange between electrons 
mediated by their interaction with a magnetic impurity 
indicates the inelastic nature of the electron scattering 
off a magnetic impurity. Indeed, using the Fermi Golden 



rule we find 



t4 pOO t4 

a{E,uj)^ J — dE'f(E')(l- f(E' +u,))oc— (2) 

^ J-oo V 



at uj T. So already in the simplest perturbation theory 
it becomes clear that there is an inelastic contribution to 
scattering. As long as E, u> T, temperature does not 
affect the inelastic cross-section in this order. It is not 
clear, however, what the relation is between the inelastic 
cross-section o~(E, uj) and the leading-logarithmic result 
Eq. QJ: On one hand, a(E,ui) oc J 4 is parametrically 
smaller than the scattering cross-section following from 
Eq. (JJJ. On the other hand, the total inelastic cross- 
section obtained from Eq. (|2J), Ctot{E) = J dus a(E,uj), 
diverges at uj — > indicating the inapplicability of the 
lowest-order perturbation results at small energy trans- 
fers. 

The lowest-order perturbation theory for K(u>,E,E') 
can be controllably improved in two respects. First, at 
E, E' ^> Tk and \oj\ <§; E, E' the four constants J enter- 
ing as a product in the perturbative result, may be re- 
placed^ by the properly renormalized 4 quantities 1 *. Sec- 
ond, the divergence at oj — ► is cut off due to the dy- 
namics of localized spin. An adequate theory may be 
developed for high temperatures, T 3> Tk, where the 
cut-off occurs due to the Korringa relaxation^. These 
improvements allow one to see that o~ t ot(E) is finite, but 
are insufficient to investigate the details of the a(E,u) 
dependence. 

In this paper we concentrate on the differential cross- 
section, a(E, uj), of inelastic scattering of a highly excited 
electron with energy E Tk- Despite the many-body 
nature of the Kondo effect, this quantity is well-defined 
at uj -C E. We show in Section [n] that in the limit 
E ^> Tk the differential cross-section is related to the 
dissipative part of the impurity spin susceptibility, x" ■ 
From the low- and high-frequency asymptotes of x" we 
extract in Sections IIIII and IIVI the behavior of the differ- 
ential cross-section o~(E,uj) in the absence and presence 
of a Zeeman energy, respectively. The analytical asymp- 
totes thus obtained are complemented by results of the 
numerical renormalization group 17 (NRG), which allows 
us to access also the intermediate range of frequencies 
and magnetic fields. The connection between the result 
of Ref. l7j and the leading-logarithmic approximation for 
the scattering amplitude (JJJ describing only elastic scat- 
tering will be explained in detail. Finally, in Section Ivl 
we discuss possible hot-electron experiments in metal- 
lic mesoscopic wires and in a semiconductor quantum- 
dot setup in order to measure the differential scattering 
cross-section ct(E,uj). 
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II. RELATION BETWEEN INELASTIC 
SCATTERING CROSS-SECTION AND 
SUSCEPTIBILITY 

The relation between the scattering cross-section of a 
"foreign" spin-carrying particle and the spin-spin corre- 
lation function of a magnetic medium is well-known from 
the theory of neutron scattering^. Here we derive a sim- 
ilar relation for scattering off a magnetic impurity of a 
high-energy electron belonging itself to the Fermi liquid 
hosting the magnetic impurity. 

The exchange interaction between the impurity spin 
and spins of electrons forming the Fermi liquid, 



k'a' 



(3) 



k.k' 



gives rise to the Kondo effect. Here J is the constant 
of exchange interaction between the impurity spin and 
itinerant electrons with energies (measured from the 
Fermi level) confined to some energy band, |ek<r| < D. 
Here s Qa / is \ times the vector of Pauli matrices. The 
Kondo problem allows for a logarithmic renormalization: 
the low-energy properties of the system described by the 
Hamiltonian © coincide with those for a Hamiltonian 
defined in a narrower band, say |eko-| < E, upon the 
proper renormalization of the exchange constant, 



J(E) 



J(D) 



1 - uJ(D) ]n(D/E)' 



J(D) = J, (4) 



where v is the density of states of itinerant electrons. 
The perturbative renormalization Eq. Q is valid as long 
as the running energy [E in the case of Eq. Q] signif- 
icantly exceeds the Kondo temperature Tjc- An impor- 
tant property of the logarithmic renormalization is that 
only exponentially wide energy intervals (61,62), such 
that vJ\ ln(£i/£2)| ~ 1 contribute significantly to the 
renormalization. That allows us to "skip" some rela- 
tively narrow strip of energies, say, (E — AE, E + AE), 
with AE <C E, in the renormalization process, yielding 
a Hamiltonian 

Hint = J(D) S ' S ««' C L c k' Q ' 

\^ a \,\e k , a ,\<D 

+ J(E) £ S-wcLtW ( 5 ) 



£-A_B<£ kQ ,e k / Q /<_E+A£; 



J(E) 



E 



S • s QQ 'C kQ c k , Q ,, 



E-AE<e ka <E+AE,\e k , a ,\<D; 
E-AE<e k , a , <E+AE,\e ka \<D 



the corresponding cross-section, within the lowest-order 
perturbation theory in J(E), can be evaluated with the 
help of the Hamiltonian ©. The first line of Eq. © 
plays the role of the Hamiltonian of a magnetic medium 
in the neutron scattering problem, and the second line 
describes the interaction of the energetic particle (we deal 
with an electron rather than with a neutron though) with 
the medium. The remaining part of the Hamiltonian 
does not contribute to the scattering cross-section in the 
lowest-order calculation. 

Consider such a scattering of an energetic electron with 
energy E and spin a in the initial and E — lu and a' , 
respectively, in the final state with oj <C E such that 
E - oj e (E - AE, E + AE). The state of the remaining 
system before and after scattering may be characterized 
by the wave functions ^ and "J f , respectively. The ini- 
tial and final state of the total system is then given by 
the product states 



initial state: \i) - 
final state: |/) 



\E, a) <g> 
\E-uu, a 1 ) 



I* 



(6) 



The differential cross-section of inelastic scattering 
a a ' cr (E, uj) is determined by the probability P a ' a (E, uj)du 
of scattering of an electron with initial energy E and spin 
a into a state within interval of energies (E — lo,E — uj — 
duj) and spin a' , 



P a > a {E,uj) duj = v F (T a < a (E,uj)dui , 



(7) 



where vf denotes the Fermi velocity. By energy conser- 
vation, uj = £/ — ^, where energies £j, £/ are associated 
respectively with the functions and '5 / involving the 
states in the domain \eka\ < D. In the absence of a mag- 
netic field, energy E is the orbital energy in the initial 
state, and to is the change in the orbital energy resulting 
from scattering. In the presence of Zeeman splitting, the 
initial energy and the energy transfer include the orbital 
and Zeeman parts, e.g. E = + og e \XBB /2. 

The standard application of the lowest-order pertur- 
bation theory in the interaction of the energetic electron 
with the remaining system yields for the scattering prob- 
ability 

Wf^t - \J(E)s aa ,(^ f \S\^ l )\ 2 2^6(i l -^ f +Lu), 

where v is the density of states for the energetic (e ~ 
E) electron. After the summation over the final states 
and proper thermal averaging over the initial states, we 
are able to relate j with a a i a (E,uj), and obtain the 
differential scattering cross-section 



with D < E — AE. The renormalized exchange constants 
here may be expressed in terms of the Kondo tempera- 
ture, vJ(e) = 1/ ln^/Tft-). There is no need to distin- 
guish between J(E - AE), J(E) or J(E + AE) as long 
as E > T K . 

If the scattering of an electron with initial energy E 
leaves it in the energy domain (E — AE, E + AE), then 



ow(£;,w) = — J 2 (E) 



(8) 



x (S a , a S zz (uj) + s+ a S+-(Gj) + s a , a S-+{u;)) , 

where sjjr v = s£, CT ± is^ l(7 . As in the theory of neutron 
scattering^, the cross-section involves a spin-spin corre- 
lation function. Here it is the correlation function of the 
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local magnetic impurity spin, 

oo 

S ab (u>)= [ dte^FW&iO)) 



(9) 



— \ /lTr_l CCShTr .\/iTr J CO 

{|*»> ,!*/>} ' 



y i \S a \y f )(y f \S b \* l )2TrS&-t f +u) 



We thus reduced the scattering cross-section to an ex- 
pression where its dependence on the energy of the scat- 
tering hot electron, E, separates from the dependence on 
the energy loss u>. The dependence on the energy loss is 
determined by the dynamics of the impurity spin charac- 
terized by the correlation function S. The spin correlator 
is related to the dissipative part of the impurity suscep- 
tibility via the fluctuation-dissipation theorem, 



(gix B fS ab {uj) = 



1 



-0U 



x'M- 



(10) 



Here /j,b is the Bohr magneton, and g is the impurity 
g-factor. The behavior of x" m various limits will be 
discussed in the following sections. The spin dynamics is 
thus included in a non-perturbative fashion. It will allow 
us to investigate the behavior of the cross-section at any 
energy transfer; at u <C Tk we apply effective Fermi 
liquid theory, and the region of intermediate energies, 
u> ~ Tk, is covered with the help of NRG calculations. 

However, it is important to note that the total scat- 
tering cross-section is fixed by the sum rule for the spin 
correlation function, S a b[oj). Consider the total cross- 
section obtained after averaging over the initial electronic 
spin configurations, tr, summing over the final ones, a', 
and integrating over the energy transfer uj, 



oo 

(Ttot(E) = J duj <7 a > a (E,id) 



3vr 1 



1 



vv F In 2 



(11) 



We substituted the explicit form for the energy depen- 
dent exchange interaction, J(E) = 1/ (v\n(E /Tk))- The 
total scattering cross-section will be used throughout the 
rest of the paper as a convenient basic unit of measure- 
ment for the differential cross-section discussed below. 

As we are mainly interested in the dependence of the 
scattering probability on the energy transfer, to, we will 
confine ourselves in the following to an analysis of the 
scattering cross-section averaged over the initial elec- 
tronic spin configurations, a, and summed over the final 
ones, a', 



a(E,u)=a tot (E) 



3tt 



S„(w) + i(5 + _(w)+S_+(w)) 



(12) 



Note that a Zeeman energy of electrons forming the Fermi 
sea was already incorporated in the definition of the en- 
ergies E and ui. The generalization of our results to spin- 
resolved scattering is straightforward. 



III. INELASTIC ELECTRON SCATTERING IN 
THE ABSENCE OF ZEEMAN SPLITTING 

In the absence of a magnetic field the expression for 
the scattering cross-section ()12f) simplifies considerably 
since the impurity spin correlator is diagonal, S(u>) = 
S zz {uj) = \S + -(uj), 



a(E,uj)=a tot (E) -S(u) 
ir 



(13) 



Let us first establish the relation between Eq. (|13|l and 
the well-known result of the leading-logarithmic approx- 
imation 3,4 ;. For that, we need to substitute in Eq. i|13|) 
the function S(u>) evaluated in the zeroth order in the 
exchange interaction J(D). In this order S^ a \uS) = 
^6(oj), which yields the well-known resuhs&i for the cross- 
section, 



r (0) 



(E,u)=a tat (E)S(u), 



(14) 



i.e., scattering is elastic in the leading- logarithmic ap- 
proximation. The elasticity breaks down, however, if one 
accounts for J(D) ^ 0. Indeed, the exchange interaction 
J(-D) leads to some dynamics of the impurity spin. The 
delta-function in Eq. Ijl4|) gets broadened, and spectral 
weight is transfered to finite energies w^O. The shape 
of the broadened peak is related to the character of the 
spin dynamics, which is different in the limits of high, 
T » T K , and low, T <C T K temperatures. We study 
the shape of the peak in these limits below. However, 
note that the broadening does not affect the total cross- 
section which is fixed by the sum rule and remains the 
same as for the elastic scattering, Eq. (|14f) . evaluated in 
the leading logarithmic approximation. 



A. Inelastic electron scattering at T > Tk 

At T ^> Tk, the local spin exhibits relaxational dy- 
namics. The Bloch equations for the impurity spin in 
the absence of a magnetic field, 



(15) 



imply the following form for the imaginary part of the 
susceptibility^, x"h( w ) = <WM with 



X M = Xo(T) — -2 

uj z + (i T K y 



(16) 



It involves the static susceptibility which is given by 
Xo(T) = (gfi B ) 2 /(4T). The decay time, t k , in the Bloch 
equations is the Korringa relaxation time 20 , 1/tk = 
ir(vJ(T)) 2 T . Inserting the scale dependent exchange in- 
teraction, J(T), the Korringa relaxation rate reads ex- 
plicitely 



1 

tk 



ttT 



(17) 
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FIG. 1: (Color online) Differential cross-section, a(E,ui), 
at large temperatures, T 3> Tk, without Zeeman splitting, 
B = 0, as given by Eq. 1191 . The Lorentz peaks have a width 
given by the Korringa relaxation rate 1/tk, Eq. 1171 . 



It is parametrically smaller than T at temperatures T 3> 
T K - 

Expression i|lb|) adequately accounts for the behavior 
of x" at low frequencies, w < T, but fails at higher fre- 
quencies. For uj ^> T, the susceptibility can be evaluated 
within the lowest-order perturbation theory in the ex- 
change constant 21 , J(D), 



i 



A i 2 M 

* win J^ 1 

i K 



(18) 



The additional logarithmic frequency dependence arises 
from the logarithmic enhancement of the exchange inter- 
action due to the perturbative RG, which is now cut-off 
at a band width D ~ u. 

The resulting differential cross-section a(E, uj) can be 
found with the help of Eq. (|10(l . It is symmetric in u> at 
small energy transfers. It shows a narrow peak at ui — 
and falls off significantly within the region of energies 
M < T: 



(19) 



where we introduced a "broadened delta-function" , 
which is a Lorentzian with linewidth T = 1 /tk , 



1 



6r(u>) = - 



1/tk 



ituj 2 + (1/t k ) 



2 ' 



(20) 



As 1/tk < T, see Eq. (|T7)). almost the full weight of the 
total cross-section is accounted for by Eqs. ljT§|) and i|2U)l. 
see Fig. ^ 

At higher energy transfers, \u>\ > T, the cross-section 
is asymmetric in to, 



a(E,uj) = a t at(E) 



l-e~"/ T to\n 2 (\u\/T K )' 



(21) 



The probability for the scattered electron to acquire en- 
ergy (uj < 0) is exponentially suppressed. Although the 
contribution of Eq. H21J1 to the total cross-section is para- 
metrically small, cx 1/ hx(T /Tk), it is worth noting that 
its decay with uj is remarkably slow. 

The slow decay of a(E,u>) vs ui is related to the de- 
pendence on the transfered energy of the cross-section 
for inelastic electron-electron scattering mediated by a 
magnetic impurity 8 . The probability for such an inelas- 
tic scattering between two electrons with initial energies 
E and E' and final energies E — uj and E' + ui was cal- 
culated in Ref. H; all of these four energies were assumed 
to be large compared to Tk- According to Ref. [see 
also Eq. (9) of Ref. 9], the contribution K (u; E, E') of a 
single magnetic impurity to this probability in the limit 
of high energy, E » |cj|, reads 



K(uj;E,E') = 
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8^1n 2 J^ 

Ik 



In 



\E'+u 



(22) 

The differential cross-section ()21|l can be obtained by in- 
tegrating K(u>; E, E') over the available phase space vol- 
ume of one of the scattering electrons, 



vf &(E, uj) 



dE'f(E')(l - f(E' + uj))K(uj- E, E'). 

(23) 



The Fermi functions in Eq. Il23[l confine the energy E' 
to an interval — uj < E' < 0. This includes a regime 
where the arguments of the E' -dependent logarithmic 
factors are not meaningful anymore and should be re- 
placed by temperature or the Korringa relaxation rate. 
After such cut-off, integration over E' is easily performed, 
yielding In -2 \uj\/Tk within logarithmic accuracy. This 
way, starting from the collision integral kernel of Ref. 0, 
one recovers Eq. (|2~T)l . 



B. Inelastic electron scattering at T < Tk 

When the temperature is below the Kondo tempera- 
ture the picture differs drastically from the zeroth or- 
der result (|14fl . For T <C Tk, the low- frequency behav- 
ior of the scattering cross-section is beyond perturbation 
theory. Nevertheless, the cross-section for small energy 
transfers, \ui\ -C Tk, may be found with the help of the 
Shiba relation 2 ^ for the susceptibility, 



(ff^ fl )VH = 27rw [xo(T = 0)Y 



(24) 



The zero-temperature static susceptibility Xo(0) is used 
conventionally 23 to define the pre-exponential factor of 
the Kondo temperature, Xo(0) = [(gfJL B ) 2 W]/ (4Tk); 
here W = 0.413... is Wilson's number. (We present 
a convenient derivation of the Shiba relation in Ap- 
pendix^]) The corrections to the Shiba relation are of 
order O (ujT 2 /T k ,uj 3 /T%) and are sub-leading. We thus 
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FIG. 2: (Color online) NRG result for a(E,cu) (solid line) 
on a logarithmic scale at T = without Zeeman splitting, 
5 = 0. A maximum at finite uj — Tk develops and scatter- 
ing with small energy transfer u) is suppressed. Whereas the 
high-frequency tail is perturbatively accessible, see Eq. 12 II 
(short-dashed line), the low- frequency tail, Eq. 11251 (long- 
dashed line), is a property of the strong coupling fixed point 
described by Nozieres' Fermi liquid theory. (W = 0.4f3... is 
Wilson's number.) The inset shows the temperature correc- 
tion according to Eq. 1251 : the contribution for negative uj is 
exponentially small for temperatures < T <C Tk- 



obtain for the cross-section at M,T <C Tk- 



a(E,oj)=a tot (E) 



W 2 



1 



to 



2 1 



1 K 



(25) 



The high-frequency limit, \ui\ 3> T%, of the scattering 
cross-section can still be obtained perturbatively and is 
given by Eq. J5TJ. 

Comparing the results Eq. (|21|l and Eq. (|25(l we see 
that for temperatures T <C Tk the differential cross- 
section cr(E, uj) peaks at energy transfers of the order 
of uj ~ Tk- It then decreases linearly upon further de- 
crease of uj, until it crosses over (at \uj\ < T) into the 
exponential tail for uj < 0, see inset of Fig. [21 At zero 
temperature the factor containing exp(— lu/T) in Eq. I|25|) 
becomes a step function which forbids any energy gain 
from the Kondo system, 

W 2 uj 
a(E,uj) = CTtot (£) e^W; (26) 
z i K 

here 9(x) = 1 if x > and if x < 0. 

The region between the asymptotes given in Eqs. I|21|) 
and H25|) can be bridged by calculations performed with 
the NRG method. In this method, after the logarith- 
mic discretization of the conduction band one maps the 
Kondo Hamiltonian onto a semi-infinite chain with the 
impurity at the end. As a consequence of the logarith- 
mic discretization, the hopping along the chain decreases 
exponentially, t n ~ A - ™/ 2 , where A > 1 is the discretiza- 
tion parameter and n is the site index. (We have used 



A = 2 throughout the calculations presented in the pa- 
per.) The separation of energy scales provided by the 
exponential decay of the hopping rate allows us to diag- 
onalize the Hamiltonian iteratively and keep the eigen- 
states with the lowest energy as most relevant ones. Since 
we know the energy eigenvalues and eigenstates, we are 
able to calculate the impurity spin correlation function 
directly (see Eq.©) given that the Dirac delta function 
appearing in the Lehman representation must be broad- 
ened when performing a numerical calculation^^. The 
result of the NRG calculation is shown in Fig. [21 

To summarize this Section we demonstrated that the 
dynamics of the impurity spin leads to inelastic electron 
scattering at all temperatures. The main contribution 
to the total scattering cross-section comes from uj ~ Tk 
or |cj| < 1/t k at T < T K and T > T K respectively. 
The total scattering cross-section is fixed by the sum rule 
for the impurity spin correlation function, see Eq. lllfl . 
and is thus determined by the effective exchange con- 
stant, J{E)i evaluated within the leading logarithmic 
approximation 3 . 



IV. ZEEMAN EFFECT IN THE ELECTRON 
SCATTERING 

We now address the case when the degeneracy of the 
impurity spin is lifted by a magnetic field. The Zeeman 
splitting of the impurity spin is described by the Hamil- 
tonian 



-gLL B S z B. 



(27) 



In the presence of the Zeeman splitting the scattering 
electron has to pay Zeeman energy in order to transfer 
spin to the Kondo system. The resonance structure for 
electron scattering involving a spin-flip will therefore dif- 
fer from the one of non-spin-flip scattering. Evaluating 
the impurity spin correlator in zeroth order in J(D) we 
obtain for the scattering cross-section (|12f> in the leading 
logarithmic approximation 



a^(E 1 uj) = a tot {E) 



1 



-0CJ 



(28) 



x - {6{uj) + 5{uj - uj z {B)) + S(uj + ojz(B))} 

The single delta- function for B = 0, Eq. (|14|> . is now split 
into three contributions. In addition to a delta-function 
at zero frequency, which is due to non-spin-flip scatter- 
ing, there are two Zeeman satellites at uj — ±ujz(B). In 
the limit of low temperatures, T <C B, the satellite at 
negative Zeeman energy corresponding to an energy gain 
of the scattering electron is exponentially small as it is 
clear from Eq. (f2"8f . 

The Zeeman energy ujz(B) depends on the renormal- 
ized g-factor, which is different from its bare value g 
appearing in the Zeeman Hamiltonian l|27|l . When we 
derived the effective interaction Hamiltonian JSJ), we in- 
tegrated out a finite band of electronic degrees of freedom 
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which lead to a renormalization of the exchange interac- 
tion J. The Zeeman term (|27H is not invariant under 
this perturbative renormalization of the Kondo model. 
Similar to the exchange interaction J, the g-factor is also 
renormalized when the band is reduced from D to D. As 
explained in Appendix[EJ the scale-dependent <?-factor in 
the leading logarithmic order is given by 



9(D) 
9 



= 11- 



1 



2 In D/T K 



(29) 



To find the observable value of g-factor, one needs to 
set D — max{T, gfisB}. The position of the Zeeman 
resonances, to the leading logarithmic order, is given by 2 ^ 



U Z {B) =g[l 



1 



-)HbB. (30) 



2 In (max{T, g^isB} /Tk 



Beyond the leading logarithmic approximation the dy- 
namics of the local spin is characterized by a further re- 
distribution of the spectral weight of the scattering cross- 
section (I28|l . However, a striking feature of the presence 
of a magnetic field is that a finite weight of the delta- 
resonance at uj = will still survive after accounting for 
the coupling of the impurity spin to the low-energy de- 
grees of freedom of the Fermi sea. In other words, at any 
ratio T/Tk a part of the scattering becomes elastic if a 
magnetic field B ^ is turned on. This can be best un- 
derstood by considering the longitudinal spin correlation 
function in time. For B ^ this correlation function will 
not fully decay with time but rather saturate at a value 
given by the finite expectation value of the impurity spin, 
(S z (t)S z {0)) -> (S z ) 2 for t -> oo. This finite saturation 
value leads to a finite weight of the delta-function 5{u>) in 
its Fourier transform and in Eq. (JSJ. Let us decompose 
a(E,u>) into the elastic and inelastic parts, 



a(E,u>) = a e \(E,u) + a inc \(E,uj) 



(31) 



The elastic part will be determined by the magnetiza- 
tion of the impurity spin 
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FIG. 3: (Color online) Weight of the elastic scattering cross- 
section, a c \(E) = J du) a e i(E, o>), see Eq. y^l ■ determined by 
NRG. The weight increases as B 2 for small magnetic fields 
and saturates logarithmically slowly to the limiting value for 
large B, see Eq. <li3"l . 



b/4t. In the developed Kondo regime, max(i, b) <C 1, the 
average spin is f(t, b) — (W/4)b. 

The weight of the elastic scattering, Eq. , evaluated 
with NRG is shown in Fig. [3J In the limit of small mag- 
netic fields this weight increases as B 2 . The saturation of 
the weight to its large-field limit, 1/3, is remarkably slow 
due to the logarithmic correction to the magnetization^, 
see Eq. ffity. 

The inelastic part of the scattering cross-section, 
(Jime^EfU)), accounts for the remaining spectral weight. 
Note however that the total scattering cross-section, i.e. 
the total spectral weight, is independent of the magnetic 
field; its value being fixed by the sum rule for the impu- 
rity spin correlator. 



a e \(E,uj) 



a tot (E)-(S z ) 2 5(co). 



(32) 



Being a thermodynamic quantity, (S z ) has a well-studied 
field and temperature dependenc o 2 ?i 2 7 . In the scaling 
regime, f(t, b) = (S z ) is a function of t = T/Tk and 
b = gfiBB/Tx- The asymptote of f(t, b) at max(£, b) ^> 1 
is with logarithmic accuracy given by 



2 ( X 21n[max(t,6)] 



(33) 



x tanh 



1 



2 ln[max(t, b)} 



Note that in the limit t = 0, b > 1, Eq. O yields the 
ground-state value of (S z ) in the perturbative regime. 
In the opposite limit of a weak field, d <C 1 « t, spin 
polarization is small according to the Curie law, / ~ 



A. Dissipative part of magnetic susceptibility 

To analyze the inelastic scattering-cross section in 
more detail for the two limiting cases T 3> Tk and 
T <C Tk, we start from presenting the proper details re- 
garding the frequency dependence of the dissipative parts 
of longitudinal and transversal impurity spin susceptibil- 
ities (x'zz an d X+-i respectively). 

At T > Tk, one may treat the exchange interaction 
J(D) perturbatively at any field B. The effect of B on 
x" is negligible as long as the Zeeman splitting ojz(B) 
is smaller than the Korringa relaxation rate 1/tk , see 
Eq. (|17|) . At higher fields, the susceptibility becomes 
anisotropic, x'zz ^ \~X-'X-i an( ^ ^ s frequency dependence 
acquires a well-resolved structure. The dissipative part 
of the susceptibility can be found from the Bloch equa- 
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tions — . The transversal part takes the form 



uj/T 2 



(34) 



where the static transversal differential susceptibility 
can be expressed with the help of Eq. I|33[) as XT = 
(gfiB)f(t,b)/B. The longitudinal part reads 



w/Ti 



where Xl 1S given by 



xl 



1 - 



ln(max{T,giJ.BB}/TK) , 



(35) 



4Tcosh 2 ^§^ 



(36) 



The factor Xl can be understood as the contribution to 
the static suceptibility which originates from the response 
of the occupation factors of the two Zeeman levels to a 
varying magnetic field, Xl — 9cfifJ-Bd(n + — n-)/dB; here 
t/cff, see Eq. (|'29p. is the appropriately renormalized g- 
factor. Note that only in the limit uoz{B) < T when the 
renormalized g-factor (|29|l is insensitive to the magnetic 
field, Xl does coincide with the full static longitudinal 
differential susceptibility xl — g^Bdf{t,b)/dB. 

In the case of a moderately high field, 1/tk <C 
luz(B) < T, the relaxation times, T\ and T 2 , equal each 
other 18 and are given by Eq. (fT7|> : T x = T 2 = t k . 

At even higher fields, loz{B) > T, the peak structure in 
X+_(w) is still described by a Lorentzian form of Eq. 134|) 
but the corresponding relaxation time is determined now 
by the Zeeman splitting rather than by temperature^, 



1 



ujz{B) 



4 In 



2 uj z 



IB) 



(37) 



The frequency dependence of the longitudinal suscepti- 
bility, however, requires additional discussion. 

Generally, the susceptibility, Xij( w )j describes the re- 
sponse of the magnetic impurity to a local magnetic field 
that oscillates with frequency to. At low frequencies, the 
variation of the dissipative part of the longitudinal com- 
ponent x" z (w) given by Eq. (|35|) can be understood in 
the framework of the Debye mechanism^ of relaxational 
losses: At co = 0, relaxation caused by the exchange inter- 
action between the local magnetic moment and itinerant 
electrons establishes equilibrium Gibbs occupation fac- 
tors for the two Zeeman-split levels. At finite but small 
frequency to, the Zeeman splitting, which is caused by the 
sum of a constant and a slowly-varying magnetic field, 
changes with time slowly, and the relaxation acts to ad- 
just the occupation factors to the instant values of the 
Zeeman splitting. The adjustment occurs via the emis- 
sion (or absorption) of particle-hole pairs with energy 
e p h ~ uiz(B) by flips of the local spin. It is the time 
variation of the occupation factors of the Zeeman-split 



levels that leads to dissipation. In the limit u> — > 0, the 
leading term in x" 2 (w), according to Eq. is 



(38) 



As was already mentioned, in the weak-field case I/Ti 
is given by Eq. ljT7|) . In the limit loz{B) > T, the time 
Ti was foundia to be T x = T 2 /2 with T 2 of Eq. gj . 
The contribution (|38|l to x" 2 from the Debye relaxational 
losses is valid at arbitrary ratio ojz(B)/T. Note however 
that despite the fact that Eq. (|38|l describes dissipation 
at low frequency the Debye mechanism is associated with 
the emission of particle-hole pairs with a comparatively 
high energy e p h ~ luz(B). In the limit ojz(B) 3> T, the 
Debye mechanism thus yields only an exponentially small 
contribution to dissipation, 



Xzz Ml Debye 



2 {g^Bfu 



ln 2 ^ 



T 



ojz{B) 



exp 



T 



(39) 

The exponential smallness of Xl come s from the small 
probability of the thermal occupation of the highly- 
excited state, corresponding to the upper of the two 
Zeeman-split levels. Temporal variations in this expo- 
nentially small quantity leads to an exponentially small 
contribution to x" 2 (w). 

Under these conditions, a second contribution, origi- 
nating from the low-energy part of the spectrum, |e| < 
max[w,T], becomes important. The processes contribut- 
ing here do not involve real impurity spin-flip processes 
(which are exponentially suppressed), but only virtual 
transitions. The starting point is the observation that 
the impurity magnetization locally polarizes the Fermi 
sea. If the Zeeman splitting of the impurity is slowly var- 
ied with a small frequency to, the magnetic polarization 
of the Fermi sea will adjust itself to the instantaneous 
adiabatic value of the impurity magnetization. Since the 
spectrum of the particle- hole pairs is continuous this ad- 
justment results in dissipation via the emission of pairs 
with small frequency e p h ~ oj, which is in contrast to the 
Debye-mechanism where the emitted particle-hole pairs 
carry a large energy of order the Zeeman splitting. As 
shown in Appendix IA1 this contribution to the suscepti- 
bility can be obtained by applying Nozieres' Fermi liquid 
theory and is adequately accounted for by the generalized 
Shiba relation Eq. l|A"6|). Evaluating d(S z )/dB with the 
help of Eq. at T <C g^BB, we find for the dissipative 
part of the longitudinal susceptibility 



7T jgn B ) u 1 
8 W |(B) ln 4 ^ 



w<w z {B). (40) 



Comparing Eq. (|4()Jl with the result for the De- 
bye mechanism, we see that the strong-field asymptote 
Eq. (|39|) for the latter mechanism is important only in a 
narrow interval of temperatures ojz(B) > T > ujz{B)/6, 
as for all practical purposes In ]xl(u)z/Tk) ~ 1. Dis- 
pensing with that interval, we will use for the dissipa- 
tive part of the longitudinal susceptibility Eq. (|35|) with 
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Ti = tk in the case u>z(B) < T, and Eq. I|40(l in the case 
cj z (B) » T. 

At low temperatures, T <C Tk, there is little effect of 
the magnetic field on x"(w) for weak fields, ff/Us-B <C Tr-. 
In the strong- field regime, u>z{B) S> 3> T, the main 
contribution to the transversal part of the dissipative sus- 
ceptibility is given by Eq. (|34l) with the relaxation time 
T2 of Eq. J37J). The longitudinal part is described by 
Eq. igni at u> <C ojz{B). Equation ||j2J) adequately de- 
scribes the non-monotonic behavior of x+_(^), but fails 
at higher frequencies; similarly the linear dependence in 
x" 2 (w) does not stretch beyond ±ujz(B). In the limit 
\uj\ ujz(B), the magnetic field does not affect signifi- 
cantly the dissipation, and Eq. lfT%|) is applicable. 



B. Elastic and inelastic components of electron 
scattering 

The coupling of the impurity spin to the low-energy de- 
grees of freedom of the Fermi seas will lead to a broaden- 
ing and redistribution of the spectral weight of the three 
delta-functions in Eq. I|28[) . 



1. High temperatures: T ^> Tk 

At high temperature, T Tk, and weak magnetic 
field, u>z{B) <C T, the spin polarization is weak, and 
the elastic component of the scattering is small. Using 
Eqs. (E2J) and we find 



0.05 



<T c \{E,Uj) = 

= a tot {E)- 



(41) 



HT/Tk) 



gfj, B B 



AT 



5(u). 



The major contribution to the scattering cross-section 
comes from the inelastic processes. At fields satisfying 
the condition uz(B)tk 3> 1, which still belongs to the 
domain of weak fields tuz(B) <C T, the single maximum 
in the w-dependence of the cross-section, see Eq. (|19H . 
splits into three: 



a inc i(E,uj) s» atot(-E) 
1 



(42) 



x - [«5 r (w) + 8 r (u - uj z {B)) + S r (uj + lo z {B))} 

The broadened delta-function was defined in Eq. i|2U|) 
with a relaxation rate, T, given by the inverse Korringa 
time, T = 1/tk- (We neglected a small part of the spec- 
tral weight which moved to the elastic component of the 
scattering cross-section). 

With the increase of the ratio u>z{B)/T, the intensity 
of the elastic scattering increases, and in the strong-field 
limit we find 



a c \{E,Lo) = <r t ot{E)^ 



1 - 



1 



ln(gfi B B /T K ) 



S(lu). (43) 




FIG. 4: (Color online) NRG result for a(E, w) at T = for 
magnetic fields g^is-B < Tk- The difference of the curves 
indicate the scattering weight which for B > is transfered 
from the inelastic to the elastic component leading to a delta- 
function peak at u) — 0, as sketched in the inset. 



Simultaneously, the maximum of o~i ne i(E, ui) at nega- 
tive u> gets suppressed, and the structure at \u>\ <C 
luz(B) broadens and becomes asymmetric. In the limit 
ljz{B)/T 1, only a single maximum at positive lo re- 
mains in the inelastic cross-section, 



a hlc i(E,uj) = otot(-E) 



3vr 1 - e-«V T w z {B) 
u/T 2 



[«-w z (B)] 2 + (l/T 2 )2 



(44) 



Here the relaxation time T 2 is defined by Eq. This 
main contribution to the inelastic scattering is propor- 
tional to X+- ( w ) an( i comes from the spin-flip processes. 
The comparison of Eqs. (gSJ and gUJ with Eq. P3|) 
shows that at luz(B) 3> Tk the effect of the dissipative 
part of longitudinal susceptibility is small starting from 
ujz(B)/T > 4. Under this condition, x"z( w ) yields a con- 
tribution to <j{E, oj) which is small compared to Eq. 144fl . 

The high-frequency tail, |w| 3> max[Tx, g/j^sB, T], is 
unaffected by the Zeeman splitting and still given by 
Eq. HU). 



2. Low temperatures: T <C Tk 

We turn now to the opposite limit of small tempera- 
ture, T <c Tk- At weak magnetic field, g^BB < Tk, the 
low-frequency behavior of the scattering cross-section is 
beyond perturbation theory. In this regime, the electron 
scatters from a fully developed, many-body Kondo sin- 
glet. Here we can use the Shiba relation Eq. 1251) to access 
the low-frequency tail of the cross-section. In the pres- 
ence of a magnetic field there are additional corrections 
to the Shiba relation of order 0(uu{gnBB) 2 /T K ) which 
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are sub-leading and are neglected in the following. We 
get for the low- frequency part \uj\ <gC Tk 



a{E,u)=a tot (E) 



W 2 



2 l-e- 

2 



j/T 



(45) 



1 / g[i B B\ u> 

6 v — J s{uj) + n 



where W is again Wilson's number—. The scattering 
cross-section decreases linearly with frequency. At uj < T 
the linear decrease crosses over into an exponential tail 
which extends to negative frequencies. In Fig.0]NRG re- 
sults at T = for the inelastic cross-section at small 
magnetic fields are compared with the NRG data at 
B = 0. In finite field the slope in the linear low- frequency 
regime is reduced. The difference in slope is of order 
0{g^sB /Tk) 2 , a correction alluded to but neglected in 
Eq. H45|) . This difference however accounts for the reduc- 
tion of the inelastic scattering weight. The weight of or- 
der 0(g/j,BB /Tk) 2 is transfered from the inelastic to the 
elastic scattering contribution leading to a delta peak at 
uj = 0, as sketched in the inset of Fig. QJ In contrast to the 
case of high temperatures (T ^> Tk^HbB) the elastic 
scattering contribution now does not sit on top of a large 
Lorentz-peak but is rather located within the scattering 
pseudogap. Although its weight is small, here it is eas- 
ily distinguishable from the background. The crossover 
from the linear dependence on uj to the high-frequency 
behavior occurs at uj ~ Tk, where the inelastic scattering 
cross-section has a maximum. The high-frequency tail is 
still given by the perturbative expression H21(l . 

When the magnetic field is increased above the Kondo 
temperature, g/J-BB 3> Tk, the elastic and inelastic 
components of the scattering cross-section are given by 
Eqs. I|43|l and l|44|l . respectively. The elastic peak at 
uj = now exhausts almost the full spectral weight of the 
longitudinal correlator, i.e. it accounts for approximately 
1/3 of the total scattering cross-section, see Fig. EI The 
remaining 2 /3 of the total spectral weight are to be found 
in the extended structure of the Zeeman satellite l|44|) 
centered at ui = ujz(B). The effect of Zeeman spitting 
on the cross-section is confined to the region of energies 
|w| < lu z {B). At \lo\ > max[T K ,gfi B B,T] the behavior 
of a(E, uj) is again given by Eq. J2U- 

In Fig. CD the inelastic cross-section is shown in the 
limit of large magnetic fields, g/igB > Tk, as given by 
Eq. g2|). The inset compares the result with the NRG. 
The low-frequency and high-frequency asymptotes are re- 
produced in the numerical calculation fairly well. The 
deviation in the width of the Zeeman peak, however, 
demonstrates the limitation of the NRG method. Due to 
the logarithmic frequency resolution the NRG tends to 
overbroaden any peak in the spectral function centered 
around a nonzero frequency. 
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FIG. 5: (Color online) Inelastic scattering cross-section for 
several magnetic field values, B S> Tk, at T = according 
to Eq. ()44|l . which is applicable unless to B. Note that the 
position of the resonance is shifted away from the Zeeman en- 
ergy as described by Eq. l.'illt . The inset shows the comparison 
of the analytical result (1441 (dotted lines) with NRG calcu- 
lation (colored lines). The low- frequency tail matches nicely. 
However, the NRG overestimates the width of the peak. 



V. POSSIBLE EXPERIMENTS 

As we have shown above, the differential scattering 
cross-section of the magnetic impurity shows a rich struc- 
ture in frequency space. In the following we suggest 
two experiments that are sensitive to the dynamics of a 
Kondo impurity and from which, in principle, the energy- 
resolved scattering cross-section can be extracted. 



A. Mesoscopic wires 

Inelastic scattering off magnetic impurities has been 
identified to be at the origin of an anomalously large 
energy relaxation in mesoscopic metallic wires^. We pro- 
pose a modification of the original experiment performed 
by Pothier et al& that allows in principle to access the 
scattering cross-section considered in this paper. We as- 
sume that the wire is connected to the reservoirs on one 
end by an open contact and on the other via a tunnel 
junction, see Fig. In the limit of small transparency 
of the tunnel junction the wire is almost in equilibrium; 
only a small amount of energetic quasi-particlcs tunnel 
into the wire and relax their energy during scattering 
processes on magnetic impurities. In lowest order in the 
transparency of the tunnel junction we can treat this re- 
laxation mechanism in terms of the differential scattering 
cross-section, a(E, uj) of a test particle coming with en- 
ergy E in an otherwise equilibrium system. 

Consider a mesoscopic wire of length L. The equilib- 
rium distribution in the right and left reservoir is given 
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FIG. 6: (Color online) Experimental setup of Pothier et aim 
with an additional tunnel barrier which limits the injection of 
hot electrons into the wire. The voltage drop across the wire 
is denoted by U. The voltage applied at the probing contact 
is V. 



by a Fermi function, f F (E) and f F (E — eU), respectively, 
where the energy E is measured with respect to the chem- 
ical potential of the right reservoir. The voltage drop 
across the wire is U . Within the wire the distribution 
function, f(E;x,U), will depend on the position across 
the wire, x € [0, L}. It is determined by the relaxation 
mechanisms and carries information on the differential 
scattering cross-section a(E,u>). This distribution func- 
tion is probed by an additional tunnel contact that is 
attached to the wire at a certain position xt € [0, L\ 
and connects it to a conductor with a sharp feature in 
the density of states. Measurement of a small tunneling 
current through this auxiliary contact as a function of 
voltage V, see Fig. [H] allows one& to probe the electron 
energy distribution. This way, the distribution function 
f(E;xT,U) in the wire at some point xt in the pres- 
ence of a bias U applied across the wire was investi- 



gate! 



]6.12 



The sharp feature in the electron density of 



states in the probe was due to its superconducting state 6 
(the BCS anomaly) or due to the Coulomb interaction in 
a low-dimensional diffusive electron system^ (zero-bias 
anomaly) . In the following we show that measurement of 
the derivative df(E; xt, U)/dU in a modified (compared 
to Ref.|a) setup of Fig. [fallows one to access the inelastic 
scattering cross-section a(E,ui). 

The distribution function within the wire is governed 
by the diffusive Boltzmann equation 29 



-D 



d 2 f(E;x,U) 
dx 2 



= I[f] 



(46) 



where D is the diffusion coefficient of the wire. The col- 
lision integral is local in space, 

/[/] = c imp v F (47) 

/OO 
dui [f(E) (1 - f(E - u>)) <j{E, to) 
-OO 

- (l-f(E))f(E-u)a(E-u,-u)] . 

where Ci mp is the impurity concentration within the wire 
and a(E,w) is the differential cross-section of a single 
magnetic impurity, and for notational convenience the 
dependence of the distribution function, /, on x and U 
has been omitted. The boundary condition at the open 



contact to the right reservoir is simply f(E; x = L,U) = 
fp(E). The boundary condition at the tunnel contact, 
which connects the wire to the left reservoir, is deter- 
mined by current conservation 

g T (f F (E-eU)- f(E;x = 0,U))= (48) 

df(E;x = 0,U) 



= -vD- 



dx 



where gr is the dimensionless conductance of the tunnel- 
ing contact and v is the density of states of the wire. 

In zeroth order in the collision integral we obtain the 
solution 

/(°>(£; x, U) = ME)^ + ME ell)-^- , 

(49) 

where we introduced the length L ; the relation of L 
to the length of the wire, L, is determined by the ratio 
of conductances of the wire and the tunneling contact, 
L /L = g w /gT, with g w = vD/L. In the limit of large 
transparency of the tunneling contact, Lq/L <C 1, the 
obtained solution reduces to the well-known formula for 
the distribution function of a diffusive wire with open 
contacts^ 9 -. However, we are focusing on the other limit 
of a large tunneling barrier, L/Lq <C 1, where we get 
x, U) = f F (E) + SfW(E; x, U) with 



5fW{E;x,U) = 

= (f F (E ~ eU) - ME)) ^ + O ( ±- 



(50) 



The deviation of the energy distribution in the wire from 
the one in the right reservoir is of first order in the small 
parameter L/Lq. 

In the following we consider the correction to the dis- 
tribution function in lowest order in the collision integral 
and in the small parameter L/Lq. We get in leading 
order in L/Lq 

I[.f {0) } - c imp v F ^ (ME ~ eU) - ME)) (51) 

/oo 
dua(E,uj) (1 - e-P") 
-OO 

(f F (E-u-eU)-ME-u>)) . 

Returning now to Eq. I|46(l . we are able to find the cor- 
rection to the distribution function. The energy depen- 
dence of df(E; x, U) /dU within the interval < E < eU 
is caused by electron energy relaxation. At T = it is 
given by 



df(E;x,U) 
d(eU) 

CimpvF / X 3 Lx 2 L 2 X 

= DLn \~6 + ~2 3~ 



(52) 



a(eU,eU-E) 
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FIG. 7: Experimental quantum dot setup. A wire leads to a 
quantum dot in the Kondo regime (indicated by the arrow), 
which is in addition weakly connected to another reservoir at 
a voltage U. The resonance-level (RL) quantum dot acts as 
a probe. 



The structure of the distribution function in this en- 
ergy interval is directly related to the differential scat- 
tering cross-section <j(E, lo). The simple relation between 
df(E; x, U) /dU and the cross-section holds as long as the 
events of scattering off magnetic impurities occur rarely 
over the time limited by the diffusion of an electron across 
the wire. Note however that in addition to Eq. 1)52 jl. 
there is a sharp contribution at the edge of the energy 
interval, E = eU, resulting from the zeroth order contri- 
bution H5()(l to the non-equilibrium distribution function. 
At finite temperature, this limits the experimental acces- 
sibility of a{E, w) for lo < T. 



B. Quantum dot in the Kondo regime 

A second experimental possibility is very similar in 
spirit to the first one but considers a quantum dot setup. 
The starting point is a semiconductor-based ballistic wire 
that has on its right hand side contact with a large reser- 
voir, see Fig. [7| On the other end, the wire is connected 
to a quantum dot in the Kondo regime. In this regime, 
the spin of the dot forms a many-body ground state with 
the electrons in the wire. Electrons injected at a bias 
U into the wire through the quantum dot form a non- 
equilibrium distribution, which is probed via an auxiliary 
weak contact having potential V. The auxiliary contact 
consists of a second quantum dot, marked RL in Fig.0 
which is tuned to the resonant tunneling regime. In the 
case of a sharp resonance, the setup of Fig. @ allows one 
to measure the electron energy distribution in the quan- 
tum wire. This non-equilibrium distribution, in turn, is 
sensitive to the inelastic transport through the Kondo dot 
and bears signatures of the differential scattering cross- 
section of the Kondo spin. 

If the left reservoir is disconnected from the Kondo dot 
the electrons within the ballistic wire have a equilibrium 
Fermi distribution, f(°>(E) = f F (E). The injection of 
hot electrons from the left reservoir will lead to a non- 
equilibrium correction to the distribution function of the 
right-movers within the wire, f(E) = f F (E)+5f(E). We 



obtain in lowest order in tunneling between the left lead 
and the dot: 

v F Sf(E) = v F J i - E) (/ F (£ - eU)(l - f F (E)) 

-e-^- E )f F (E)(l-f F (t-eU))), 

(53) 

where we used already the detailed balance relation, 
a{E,Lo)e"^ w — a(E~uj, — lv). After taking the derivative 
with respect to U , the above equation simplifies consid- 
erably at T = 0, and we get 



dm 

d(eU) 



a{eU, eU - E) . 



(54) 



The measurement of this quantity with help of the aux- 
iliary contact thus yields direct access to the differential 
inelastic scattering cross-section of a Kondo system. 



VI. SUMMARY 

We analyzed inelastic scattering of energetic electrons 
off a magnetic impurity. For such scattering, the depen- 
dence of the differential cross-section, o(E, lo), on energy 
E of the incoming electron is logarithmically weak at 
E 3> Tk, and arises from the renormalization of the ex- 
change coupling. In the leading logarithmic approxima- 
tion, the total cross-section a tot — J dujcr(E, lo) is propor- 
tional to 1/ Iv?(E/Tk), in agreement with Ref. 0. More 
interestingly, the electron scattering is inelastic, and the 
dependence of <j(E,oj) on the energy transfer, lo, is de- 
termined by the spin-spin correlation function of the im- 
purity or, equivalently, by the dissipative part of the im- 
purity spin susceptibility, \" ■ In the absence of magnetic 
field, the elastic component of scattering appears only in 
order 1/ In 4 (E/T K ). 

Our findings confirm and quantify the conclusion of 
Ref. IjJ regarding the inelastic nature of Kondo scattering, 
and also provide a clear physical picture of the mecha- 
nism of inelastic scattering. In the absence of magnetic 
field, the inelastic scattering cross-section is parametri- 
cally larger than the elastic one. The typical energy 
transfer \lo\ in an inelastic scattering event is however 
small compared to E. At high temperatures, T » Tk, 
the characteristic energy transfer is determined by the 
Korringa relaxation rate of the magnetic impurity, and 
at low temperatures, it is defined by the value of Tk- In 
the high-temperature limit, the cross-section is maximal 
at lo — 0, and at T <C Tk it reaches its maximum at 
lo ~ Tk ■ The decrease of the cross-section in the domain 
lo Tk is remarkably slow, a(E,u>) cx [lo In 2 (lo/Tk)]^ 1 ■ 
The domain of intermediate energy transfers, to ~ Tk, 
is covered by NRG calculations. The numerical results 
fit well with the analytically evaluated asymptotes at 
lo <C Tk and lo 3> Tk- In the presence of an exter- 
nal magnetic field, the Zeeman splitting of the magnetic 
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impurity levels results in the appearance of an elastic 
component of electron scattering already in the leading 
logarithmic order (in E/Tk), 

Finally, we proposed possible hot-electron experiments 
with a metallic mesoscopic wire and with a semiconduc- 
tor quantum-dot device which in principle allow one to 
access the differential scattering cross-section of a local- 
ized magnetic moment. 
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APPENDIX A: DERIVATION OF SHIBA 
RELATION 

Here we provide a simple derivation of Shiba relation^, 
using Nozieres' idea of a low-temperature Fermi liquid 
description of the Kondo problem. 

Within Nozieres' theory, at T <C Tk the effect of a 
weak [g[iBB -C Tk) magnetic field applied to a Kondo 
impurity is described by a local-field Hamiltonian, 



Hr 



2\ 



gu. B v 



(Al) 



k,k' ,cr 



Here v is the density of states at the Fermi level, xo — 
[W / (5/is) 2 ]/(4Tif) is the linear susceptibility, summation 
over fc and fc' occurs within a shell of states Afc suffi- 
ciently close to the Fermi level (Afc ~ Tk/vf with vf 
being the Fermi velocity), and field B is applied along 
the z-axis. One may easily check that the action of the 
field described by the Hamiltonian (|A1(I indeed results in 
a local magnetization M = XoB- For that one starts with 
the evaluation of the spin-dependent scattering phase S a 
off the local perturbation Eq. IjAlfl using the Born ap- 
proximation, 



S a = 



Xo 
9VBV 



B. 



a = ±1. 



(A2) 



Having the phase difference 8^ 



we evaluate the 



magnetization using the Friedel sum rule, 
gUB 5+ - 5- 



M = 



XoB. 



(A3) 



Having the right form of the local perturbation, we now 
allow for a slow variation of the field, B = Bo cos(wi), as- 
suming that the frequency u> <C Tk- Next we evaluate the 




FIG. 8: (Color online) NRG comparison of the low-frequency 
asymptote of x'zz with the prediction of the generalized Shiba 
relation Eq. ijABfl (dotted lines). The inset shows the nu- 
merically evaluated magnetization whose derivative enters 
Eq. 



energy absorption rate w caused by such time-dependent 
perturbation. Using the Fermi Golden rule, we arrive at 





XoBq' 




W = TTUJ 


gu, B v _ 





v 1 / def(e)[f(e-uj)-f(e + w)] 



XoB 



9HB 



(A4) 



In the last line, we discarded corrections of order 
0(eT TK > T ) arising from the boundaries of the energy in- 
tegral. Recalling finally that w — \ux" {u)Bq, we arrive 
at the Shiba relation Eq. . 

Using the framework of the above derivation, it is 
straightforward to generalize the Shiba relation to the 
case of a weak slowly varying field applied to the local 
moment on top of a time-independent field B of arbitrary 
strength. In the generalized relation, xo is the static dif- 
ferential susceptibility, and the relation is applicable in 
the regime uj,T <^ max{_B, Tk}- 

We assume that the basis of the effective low-energy 
Hamiltonian has been chosen such that it incorporates al- 
ready the effect of the time-independent local magnetic 
field B. Consider now a small perturbation to this ef- 
fective Hamiltonian induced by a small change in the 
applied local magnetic field B + SB, 



H SB = 



2 DM 
gi^BV dB 



SB 



^ S aa c \<j C k' (? ■ 



(A5) 



The summation over fc and kl is bounded by |fe|, \k'\ < 
max{t;/iB_B, Tk}/vf- The prefactor can be determined 
in the same way as before. In contrast to the limit 
B = 0, here the resulting phase shift yields information 
about the change in magnetization M(B + SB) — M(B) = 
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(a) -vW^Ov (b) ( c ) rS^\^J^\ 

FIG. 9: (a) First order Knight-shift diagram; (b) two- 
loop self-energy correction; (c) one-loop vertex correction. 
The wiggled line represents the propagator of the Abrikosov 
fermions and the solid line the electron propagator. The dot 
signifies the Kondo interaction. 



(dM/dB)5B, where dM/dB is the differential suscepti- 
bility. The same arguments as above will yield the gen- 
eralized Shiba relation 



d(S z 



dB 



(A6) 



Here (S z ) is the equilibrium average spin value in the 
presence of field B. In the perturbative regime, the av- 
erage is given in Eq. I|33l) . 

In Fig. [S] the prediction of the generalized Shiba rela- 
tion (|A6|) is illustrated with the NRG result. The sus- 
ceptibility, Xzzi an d the non-linear static susceptibility, 
Xo = 9^Bd(S z )/dB, have been independently evaluated 
with the NRG. The dashed line in Fig. is plotted with 
the help of Eq. (| A6|) and compares well with the low- 
frequency asymptote of \"z ■ 



APPENDIX B: RG EQUATION FOR THE 
IMPURITY g-FACTOR 

We present a derivation of the two-loop RG equation 
for the impurity (/-factor of the Kondo model and, in 
particular, explain the different roles played by the im- 
purity and conduction electron (/-factor in the renormal- 
ization process. To this end, we will use Abrikosov's 
pseudo-fermion representation^ for the impurity spin, 
S = /tier/, where p = (/|,/|) in a compact spinor 
notation and a is the vector of Pauli matrices. We will 
need the action of the Kondo model, which consists of 
three parts, S = S s +Sd + Sx- The quadratic part of the 
Abrikosov pseudo-fermions reads 



S d 



P 



dr/t(r) 



d T - A - g^ B Tj^Ba 



f(r), (Bl) 



where g is the impurity g-factor and B a is the magnetic 
field, which is taken to point in the z-direction B a = 
B8 az . In order to enforce the Hilbert space constraint, 
f< f = 1, a chemical potential, Ao — * oo, is introduced 3 . 
The Kondo interaction is given by 

Sk = J**r U\t)\<t^{t)\ J ab (V(T)i<r 6 /(T)) , 

(B2) 

where the local electron operator at the impurity 

site is *t = / 2^( c !t' c !i)- We allow for dirfer - 
ent values of the exchange interaction in the direc- 
tion orthogonal and perpendicular to the magnetic field, 



(Jab) =diag{ J±, Jj_, Jjj}. Finally, the quadratic part of 
the s-electrons reads 



D 



S s =f Q dr J 2^4>) [fir + k] c ka (r) . (B3) 

-D 

In the presence of a Zeeman energy for the s-electrons, 
the Fermi sea of the spin-up and -down electrons are 
shifted with respect to each other giving rise to a finite 
Pauli magnetization. In Eq. (|B3|1 we assumed that the 
band has already been symmetrized around the respec- 
tive Fermi energies by integrating out a finite number of 
electronic degrees of freedom. This process results in a 
perturbative renormalization of the impurity g-factor due 
to the so-called Knight shift. The first-order Knight-shift 
diagram is shown in Fig. [5^. The (/-factor, g, appearing 
in (|_B1I) is therefore understood to be already the Knight- 
shifted impurity g- factor, 



J\\v 



9 = 9i 



0{J\\vf, 



(B4) 



where gi and g e are the bare impurity- and electronic 
g-factors, respectively, and the density of states is v — 
1/(27ti>f). As is clear from Eq. (|B4|> the electronic g e 
can be absorbed in an effective impurity (/-factor. The 
Knight shift is thus only a perturbative phenomenon and, 
in particular, is not enhanced by logarithmic renormal- 
izations. This is expected since the Pauli magnetization 
affects only electronic states far away from the Fermi edge 
deep inside the Fermi sea. 

The field theory can be renormalized-Q with a wave- 
function, impurity (/-factor and Kondo-coupling renor- 
malization (in addition to a counterterm absorbing a shift 
in the unphysical chemical potential Ao), 



f = ^Zf R , g = 



Jab = 



jR 

J ab 



(B5) 



We compute the renormalization of the Kondo coupling 
to one-loop and the renormalization of the wave-function 
and (/-factor to two-loop order. The corresponding di- 
agrams are shown in Fig. ^jp and c. The resulting RG 
equations for the Kondo vertex are the well-known poor 
man's scaling equations 4 



d(J±v) 
d\nD 



-(J±v)(J\\v) , 



d(J\\v) 
d\nD 



-{J xv) 



For the wave-function renormalization we obtain, 
d\nZ 1 



d\nD 8 



((J^f + 2(J^f) 



(B6) 



(B7) 



Finally, the main result is the RG equation for the g- 
factor, 



d 9 1 / j x2 

= t;9(J-lv) . 



dlnD 2' 



(B8) 



15 



Solving this equation in the isotropic case, Jj_ = J|| = J, 
and expanding the result in leading logarithmic order we 
get 

■'<°» =4- ^* + *) t) <b9 » 

where we already substituted the Knight-shifted g-factor 
l)B4[l. In the scaling limit, J — > while Tk is held fixed, 
any dependence on the electronic ^-factor, g e , vanishes, 



and we obtain the result cited in the body of the pa- 
per, Eq. (|29H . In particular, note that in the absence 
of a Knight shift, g e = 0, the perturbative correction to 
the ^-factor starts only in second order in the exchange 
coupling J but, nevertheless, after renormalization group 
improvement leads to a correction which is of leading log- 
arithmic order. At zero temperature the RG equation for 
the g-factor coincides with the RG equation for the im- 
purity magnetization, which was already determined in 
Ref. EH 
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